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ABSTRACT 


The heavy-lift launch vehicles (HLLV) proposed for use in constructing satellite power systems (SPS) 
would deposit various contaminants in the middle atmosphere, contaminants that could conceivably have 
adverse effects on climate and upper air structure. These contaminants consist of the major constitu- 
ents of water vapor, hydrogen, carbon dioxide, and carbon monoxide, and the minor constituents of sulfur 
dioxide and nitric oxide in the rocket effluent, as well as nitric oxide formed during reentry. 

To assess the magnitudes of the effects, we have constructed new models or modified existing models. 

They are: one- and two-dimensional photochemical models, a plume model, a noctilucent cloud and contrail 

model, a reentry model, and a model of the lower ionosphere, all of which are described in detail in the 
report. Using a scenario of 400 launches per year for 10 years, we have performed assessments and 
arrived at the following conclusions: (1) The buildup of water vapor, nitric oxide, CO2 > and CO, and 

sulfur dioxide, including a possible "corridor" effect (zonal enhancement centered on the launch lati- 
tude) is not likely to be significant, (2) Perturbations to odd-oxygen (0, O3) are not likely to be sig- 
nificant in that ozone total column density decreases are probably less than 0.1%; in fact, ozone 
increases due to NO deposition may occur. (3) Perturbations to odd -hydrogen species (H, OH, HO2) are 
not likely to be significant in the stratosphere and mesosphere; however, there could be a doubling of 
thermospheric hydrogen. (4) The effect of water vapor emission on global cloud coverage would be negli- 
gible. Although the possibility of impressive local noctilucent cloud displays near the launch latitude 
is possible, we do not expect them to have any climatic significance. (5) Perturbations of HF, VLF, and 
ELF communication links will be minimal compared to normal variability, however, D-region electron con- 
centrations and conductivity may be greatly enhanced locally during atmospheric reentry. 
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1 . INTRODUCTION 


The U.S. Department of Energy is responsible for examining new or alternative energy sources. One such • 
potential energy source which has been proposed is that of solar power collected by suitable satellite 
systems and transmitted to Earth in the form of many microwave beams. Among the important factors, 
which must be addressed as part of the feasibility investigation for such a system, is the environmental 
impact of the heavy-lift launch vehicles (HLLV) necessary for the construction of solar power systems in 
space. The purpose of the study reported here is to assess the effects on the middle atmosphere (strato- 
sphere, mesosphere, and lower ionosphere) of launching about 400 HLLV*s per year for 10 years or more; 
other investigating teams are dealing with environmental effects on the troposphere, thermosphere, and 
upper ionosphere. 


Table 1. HLLV Launch Effluents (Ref. 1) 


The principal emissions from the HLLV that can be expected to modify atmospheric composition are water 
vapor, hydrogen, and carbon dioxide deposited by rocket launch engines (Table 1), and nitric oxide formed 
during reentry. In addition, small amounts of nitric 
oxide and sulfur dioxide are expected to be formed in 
the combustion process and deposited during launch. 

Water vapor and hydrogen deposition profiles given in 
units of molecules cm”^ are presented in Fig. 1; note 
the discontinuity at an altitude of 56 km, where the 
first stage, which uses methane as fuel, cuts off and 
the second stage, which uses hydrogen as fuel, 
ignites. Also note the large Increase in the deposi- 
tion rate as the vehicle approaches an altitude of 
120 km; the increase is a result of the nearly hori- 
zontal flight near the end of the trajectory and 


Effluent 

Mass per launch, 
tonnes 

Altitude interval, 
km 

H 2 O 

3860 

10-125 

H 2 

160 

10-125 

CO 2 

1650 

10-56 

CO 

750 

10-56 

SO 2 

<0.5 

10-56 


burnout at 120 km. Carbon dioxide, a product of 
methane combustion, is also emitted up to first 
stage burnout; its deposition rate has been 
estimated to be 0.4 of the water vapor rate 
(Ref. 1). The emission rate of nitric oxide 
during the launch phase is also shown in Fig. 1. 
The emission rate of sulfur dioxide is estimated 
by assuming that the first-stage fuel contains 
0.05% sulfur.^ Because the production of nitric 
oxide during reentry is the subject of a rather 
lengthy calculation, we defer its discussion to 
a later section. 

Because of the complexity of the physical pro- 
cesses occurring in the stratosphere, mathe- 
matical models are required to obtain quantita- 
tive estimates of the effects of HLLV operations. 
General models that consider the interactions of 
radiation, chemistry, and dynamics are beyond 
current capabilities. Hence, a range of more 
specific models must be used to describe each 
area of concern. 



H2» H2O DEPOSITION, molecules cm 
NO DEPOSITION, molecules cm 


-1 


X 10 


23 


1 X 10^0 


Fig. 1. 


Deposition rates (molecules cm”^) of water 
vapor, molecular oxygen, and nitric oxide assuming a 
launch rate of 400 heavy lift launch vehicles per 
year. Note the staging at 56 km. 


In the following sections we discuss the atmo- 
spheric photochemical models used in this 

assessment. Including a model for high-altitude cloud formation, a calculation of nitric oxide produc- 
tion during reentry, and an ionospheric model. We then apply the models to the assessment of HLLV 
effects on the middle atmosphere, evaluating both possible global effects as well as short-term local 
effects, and we discuss uncertainties in the assessment. Finally, we summarize our conclusions, Identi 
fying those effects that could be significant and may require further study. 


2. ATMOSPHERIC PHOTOCHEMICAL MODELS 


The atmosphere is extremely difficult to model. Even the tenuous stratosphere and mesosphere present 
great difficulties. In fact, atmospheric simulation by mathematical models has been called the "last 


^R. L. Brubaker, Argonne National Laboratory, private communication, 1980. 
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great unsolved problem of classical physics." Such statements usually refer to the simulation of atmo- 
spheric motions, but even the considerably less complex problem of simulating atmospheric photochemistry 
introduces a number of approximations that depart from reality. For example, the set of chemical species 
and reactions that can be accommodated is necessarily restricted in scope, atmospheric variability cannot 
be simulated, and large-scale long-period averaging must be invoked. 

The transport of species is simplified to an even greater degree. One-dimensional models, which in some 
sense simulate average global conditions, can represent only vertical transport and usually do so with 
the so-called "eddy dif fusivity. " Two-dimensional models are less limited in that they can include an 
"eddy-diffusion tensor," which represents motions leading to mixing, and meridional bulk motions. In 
the model used in the study reported here, the bulk velocities are obtained by zonally averaging the wind 
fields predicted by three-dimensional models; eddy coefficients were obtained by an adjustment procedure 
based on matching computed tracer distributions with their observed distributions. 

To carry out the required assessments most efficiently and realistically, two photochemical models are 
used: a one-dimensional (1-D) model and a two-dimensional (2-D) model. The 1-D model, which contains 

over 50 species, is most useful for investigating globally or hemispherically averaged changes in compo- 
sition; it is also capable of simulating a much larger system of species and reactions than the 2-D 
model, which is limited to about 25 species below 55 km and 9 species above 55 km. The 2-D model, oh 
the other hand, can simulate meridional variations in composition changes and is thus able to treat pos- 
sible seasonal and "corridor" effects; the latter refers to enhanced changes in a zone containing the 
launch or reentry windows. Because both models have been described in detail in the literature 

2—5), only brief descriptions are given here. Both models are based on the continuity equations 
for trace constituents 


3t 




( 1 ) 


where and are production and loss rates, respectively, for the ith constituent, v is the 
meridional bulk velocity, is the "eddy flux" representing large-scale eddy motions, and n^ is the 
number density of the ith constituent. In the 1-D model the velocity v is, of course, absent, and 
the eddy flux is written as 



( 2 ) 


where e^ is a vertical unit vector, is the eddy dif fusivity, and T and H are the temperature and 
mean atmospheric scale height, respectively. Because the 1-D model extends into the lower thermosphere, 
molecular diffusion cannot be neglected; a term representing molecular diffusive flux must there- 

fore be added to 



Here, is the scale height of the ith constituent 


r kT ra 

[ = = — 

i mg m 


(3) 


(4) 


mi and m are the molecular mass of the ith constituent and mean molecular mass, respectively, k is 
the Boltzmann constant, and g is the gravitational acceleration. The 2-D model, which does not extend 
into the thermosphere, ignores molecular diffusion. However, the eddy diffusion is more complicated 
than in the 1-D model because it is a tensor, that is, it contains "off-diagonal" elements, 

= + . (5) 

which are Important below 25 km. 
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2.1 The 1-D Model 


The 1-D model extends from an altitude of 10 km to 120 km with a 2-km grid spacing; time steps are vari- 
able, depending on the nearness to equilibrium. The form of the flux equation used in the model is 
slightly different from that given above in that mixing fractions rather than number densities are used. 
The flux equation, which now contains both eddy and molecular diffusion terms, is 


where is the mixing fraction. 



( 6 ) 


(7) 


and Tq and Mq are atmospheric temperature and number density, respectively, at a reference altitude, 
Zq. The vertical profile of eddy diffusion used in the model is shown in Fig. 2. 


Solutions of the species continuity equations, 
which are based on approximate, but very stable 
and efficient, numerical integration schemes, 
often violate mass conservation within groups 
of related species such as the nitrogen oxides 
(NO, NO 2 , NO 3, N 2 O 5 , HNO3). This problem is 
overcome in the model by employing a technique 
that accurately monitors the total concentra- 
tions of nearly conserved sets, or families, of 
compounds (Ref. 2); the concept of aeronomic 
families is well-known in aeronomy. Our tech- 
nique uses an equivalent continuity equation 
for an entire family of gases, which is obtained 
by summing the continuity equations of all the 
family members after each is multiplied by an 
appropriate factor, a^, representing the number 
of odd-atoms of the type considered in the 
species. The criterion for assigning a is 
based on the time constant for switching odd- 
atoms among species; if the time constant is 
short a > 1; but if it is long, a = 0. Thus, 
summing ( 1 ) for a particular set of species we 
obtain: 


8t “ ^SJl ” ^SJl ” 3z ’ *■ = 1. 2, . . . L , (8) 

where i is the family index. 


P 


SI 


S|, " 2 . (9) 

li), 


53 - P.L.) . 

U}, 

(10) 


(11) 



1q3 10^ 10^ 

EDDY DIFFUSION COEFFICIENT, cm^ sec"'^ 


and {i}^^ indicates the subset of species in the Fig. 2. The eddy diffusion coefficient (K) profile 
S-th family. used in our one-dimensional model studies. 
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The family flux, which follows from (6), is 


'H- y 3z Hy \m 7 

The factor mj^/m in (12) is a weighted average value. 


(K + D) 

y 


( 12 ) 



Xi “i^i (^)/X “I'l ’ 

U}„ {i}„ 


(13) 


which is nearly time invariant, because the ratios of species concentrations within a family are nearly 
constant under most conditions. Obviously, the species boundary conditions, which we discuss later, are 
also additive for families. 


What makes the family technique useful is that the net photochemistry for the family, represented by 

(S) > usually contains only slow processes that affect the total family abundance, and pre- 
cludes the rapid cyclic chemical interactions that occur between family members. Accordingly, the summed 
continuity equations in (8) are inherently more stable than the individual species continuity equations, 
and their solution can be made very accurate. 


Family concentrations obtained from (8) are used to correct the species abundances computed using (1) 
and (6). The families of species utilized in our model are given in Table 2 in the hierarchical order 
in which they are solved and applied to correct the species concentrations. For example, HNO 3 is a mem- 
ber of the hydrogen and nitrogen families, but it is adjusted during the hydrogen correction cycle and 
remains fixed during the nitrogen cycle. 


Hydrogen: 

Carbon: 

Nitrogen: 

Sulfur: 

Chlorine-oxide : 
Odd-chlorine : 
Bromine: 

Oxygen: 


Table 2.- Families of Species 


H(l), 0H(1), HOad), H202 (2), HN02(D, HNOsd), HOCl(l)^ 

CH(1), CH2(1), CH3(1), CHO(l), CH20(1), CHaOCl) , CH302(1), CH402d) 

N(l), N0(1), N02(1), NOad), HN02d), HNOad), N20 s( 2), C10N02(1)^, Br0N02d)^ 
S(l), S0(1), S02(l). SOsd), HS03<1), H2S04(1) 

Cl(l), CiO(l), C102d), ClONOzd), HOCl(l) 

Cl(l), ClO(l), C102d). ClONOod), HCl(l), HOCl(l) 

Br(l), BrO(l), HBr(l) , BrONOaCl) 

0(1), 03 d), O^D(l), NOad), N0s(2), HNOad), N 20 s( 3 ), ClO(l), C102(2), C10N0?(2), 
BrO(l), Br0N02(2) 


^The number in parentheses is the species weight a within the family. 

The species is not corrected in this family, but rather in the corresponding halogen family. 


Two types of families are easily recognizable: one in which the family members recombine in pairs 

(hydrogen, nitrogen, oxygen), the other in which the members recombine with species outside of the fam- 
not recombine at all (the other families). Families can also be classified according to 
their internal chemistry: the members of some families cycle rapidly among themselves in reaction loops 

(hydrogen, nitrogen, oxygen, chlorine), while the members of other families react in a chain, one species 
leading to the next, from the family source to its sink (carbon, sulfur). Note that the consideration of 
families requires the identification of the atom carriers in key aeronomic processes. For example, in 
the oxygen family we can readily find reactions that might not seem to involve oxygen atom transfer or 
recombination, but actually do; thus, the reaction of CIO with NO to form Cl and NO 2 (which is rapidly 
photo lyzed to NO and 0) does not involve any odd-oxygen production or loss according to our criteria for 
odd-oxygen (see Table 2). Thus, the family concept is helpful in understanding aeronomic reaction 
mechanisms. 


The chemical rate coefficients used in the model are essentially those tabulated in Ref. 6, except that 
the photolysis rate for NO is taken from Refs. 7 and 8. Furthermore, we include a source of mesospheric 
N( D) (first excited state of N) due to the ionization of NO by solar Lyman-a radiation: 

NO + hv^ No"^ + e , (14) 
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(14) concluded 


NO*^ + e -> N(^D) + 0 . 

In the model, diurnal averaging of photodissociation rates and reaction rate coefficients is performed 
using, respectively, the techniques reported in Refs. 9 and 10. 


For each species we specify a lower boundary flux and an upper boundary flux, which may be 

fixed or have a specified time dependence. At the lower boundary we also include a flux component. 




(15) 


by defining the "velocity" at the boundary, vg^, and concentration, ng^ (if ng^ = 0 , our code sets 
“ 0 automatically); nj^^ is computed by the model. Whenever we use boundary condition (15), we 
usually set the boundary velocity, vg^, equal to 1 cm sec"^, which is estimated from typical tropospheric 

mixing rates. A justification for this approach together with its mathematical application is given in 
Ref. 3. Upper boundary conditions are imposed in a somewhat similar manner except that allowance must 
be made for the escape of hydrogen. The mean escape flux of H atoms is about 10^ cm“^ sec"^ (Ref. 11), 
which must be included in the computations of both H and Hz. We specify the upper boundary fluxes as 


= bnjj - , | 

I (16) 

% = ) 

with the total escape flux ' given by 

The constant b is then determined from simul- 
taneous estimates of the escape flux and nj^, 
while a must be estimated from the flux of Hz 
at 120 km. 

In Fig. 3, 1-D model predictions for the ambient 
atmosphere are compared with observations. The 
agreement between prediction of atomic oxygen, 

0( P) , and its observation is reasonable consid- 
ering the large zenith angle of the observation; 
the high-altitude computed concentrations are 
perhaps a factor of 2 too large between 90 and 
100 km, but at other altitudes the agreement is 
quite good. Computed ozone concentrations are 
generally within the error bar of the Krueger- 
Minzner (Ref. 16) empirical model, although the 
high-altitude (above 45 km) predicted concentra- 
tions are somewhat too small.- Predictions for 
OH and water vapor appear to lie within the range 
of observed values; however, it should be remarked 
that the range of observed water vapor is quite 
broad, extending from as low as 2 ppbv to 10 ppbv. 
Also, the OH observations were made at large solar 
zenith angles whereas the computed values are for 
average mid-latitude daytime conditions; to cor- 
respond to the same conditions, the predicted 
values would have to be reduced by a factor of ~ 2 . 
The predicted nitric oxide abundance in the 
stratosphere lies well within the range of mea- 
surements. However, in the upper mesosphere and 
lower thermosphere, there is disagreement; this 
difference may be due in part to the large zenith 
angle of the observations and in part to adjust- 
ments needed in the mesospheric NO photochemical 
sink and in the diffusion rate. 



Fig. 3. Vertical profiles of H, OH, HzO, 0(^P), O 3 , 
and NO calculated with the one-dimensional model. 

The concentrations are average daytime values. Mea- 
sured concentrations shown are; 

OH: \-m~\ Anderson (Ref. 12); White Sands, N. Mex. , 

in April with solar zenith angle of 86 ® 
l-i^l Anderson (Ref. 13); Palestine, Tex., in 
July and January with solar zenith angle 
of 80® 

HzO: 0-0 Summary and review of measurements given 

by Ellsaesser et al. (Ref. 14) 

0(^P); (~^I Hewlett et al. (Ref. 15); White Sands, 

N. Mex. , in December with solar zenith 
angle close to 90® 

l-o-i Summary and review of measurements given 
in Ref. 6 

O 3 ; l-Ar*[ Krueger and Minzner (Ref. 16) 

NO: x-x Summary and review given in Hudson and 

Reed (Ref. 6 ) 

T Baker et al. (Ref. 17); White Sands, 

N. Mex. , in October with solar zenith 
angle of 92® 
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2.2 The 2-D Model 


The 2-D model extends from the surface to an altitude of 90 km in 2.5-km steps and from lat. 80°N to 
lat , 80°S in 5° steps; time steps are fixed at 1 day. The form of the flux equation used in the model is 
that given in (5) above. The temperature T, which is a function of season and latitude as well as 
altitude, is taken from Ref. 18 for altitudes below 20 km and from Ref. 19 for higher altitudes. The 
two-dimensional divergence operator is represented in the meridional plane of a sphere by 


V« 


iR cos e 90 ^ 



where 0 is the latitude, e 0 is the tangential unit vector, and R is the radius of Earth. 


(18) 


Transport parameters are specified in the model for each season. They are of two types: bulk velocity, 

which represents large-scale mean meridional motion, and "eddy dif fusivity which represents smaller 
scale motions. The parameters that we have adopted are physically reasonable and lead to satisfactory 
predictions of observed tracer distributions. 


The mean meridional circulation is obtained by the kinematic method from the averaged equation of mass 
continuity, which states that 


•“ + V • (pv) = 0 , (19) 

where p is the atmospheric bulk density and v is the large-scale atmospheric velocity. With the 
assumption that the density field is in a steady state, the approximate form of this equation in spheri- 
cal coordinates is 


1 3 . 3(pw) 

^ (pv cos e) + = 0 


( 20 ) 


where the overbar denotes an average with respect to time and longitude, and v and w are the mean 
meridional and vertical velocity components. Equation (20) implies the existence of a "stream function" 
for the total flux such that 


2 i 7 Rpv cos 0 = - 1^ ; 27 tRpw cos 6=4-^* (21) 

o Z Rod 

If the distributions of p and v are known, the first of equation (21) can be integrated vertically to 
obtain and w can then be obtained from the second. 

The density was obtained from the hydrostatic equation using temperatures discussed previously. The 
mean meridional circulation is based on that published by Cunnold et al. (Ref. 20), which was extrapo- 
lated upward to 90 km. The resulting stream function satisfies (21). 

The eddy fluxes are modeled by diffusion coefficients Kyy> Ky 2 , and K 22 (the local "y" axis is tangent 
to the meridian; i.e., dy = Rd) . The magnitude of K 2 Z this model is generally somewhat smaller 

than that in a one-dimensional model because vertical transport in two-dimensional models is affected by 
large-scale bulk motion (w) as well as by eddy mixing. The assumed values of Kyy, which is the horizon- 
tal analog of K 22 > are strongly latitude-dependent, but height-independent. 

The most interesting of the j , however, is the "off-diagonal" Ky 2 (=K 2 y) because it leads to an 
apparent counter-gradient flux of ozone. It has been convincingly argued by Newall (Ref. 27) that the 
"mixing surfaces" in the lower stratosphere slope more steeply toward the pole than do the surfaces of 
potential temperature. That is, parcels of ozone moving poleward in the lower stratosphere tend also to 
move farther downward in the transport region to altitudes where their destruction becomes less effi- 
cient. Hence we see a buildup of ozone column density at high latitudes. Such an effect is associated 
with diabatic processes, which cause the relative slopes of surfaces of constant mixing and constant 
potential temperature to shift. One can in principle calculate the by computing the variances of 

observed winds over a long period of time. In practice, one can obtain only order- of -magnitude estimates 
in this way. Our , which are based on observed tracer distributions, are consistent with the order- 

of-magnitude estimates obtained from wind-variance calculations. In order to facilitate computation of 
the j , we employ a suggestion by Brasseur (Ref. 22) that they be expressed in the form 
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( 22 ) 


''ij = • 


where 0 is a function of latitude alone and Z is a function of height alone (note that our Zyy«l), 
These functions are readily adjusted to yield K-fields that result in smooth distributions of trace 
constituents. 

End boundaries are taken at lat. 80 ®S and lat. 80 °N because meridional fluxes are expected to be small at 

those latitudes. Hence, the end boundary conditions are taken to be zero flux of all constituents 

across the vertical boundaries. The upper boundary conditions are given by setting the fluxes equal to 
zero for all species except H, NO, and 0(^P). We have checked the upper boundary conditions (at 90 km) 
for the various species against results from a one-dimensional model that extended up to 120 km. It was 
found that the effect of choosing a mixing equilibrium condition at the upper boundary had very little 
effect on any of the constituents at altitudes below 50 km. The lower boundary condition used for all 
species except N 2 O, CH 4 , HNO 3 , NO 2 , O 3 , HCl, H 2 O 2 , and the halocarbons, is chemical equilibrium because 
of their short lifetimes against chemical loss. Because HNO 3 , HCl, and H 2 O 2 are water soluble, their 

number densities are set equal to zero at the lower boundary. The number densities of CH 4 , N 02 > and N 2 O 

are fixed at 3.7 x 10^^ cm“^, 3 x 10^ cm“^, and 7.5 x 10 ^^ cm“^, respectively, at the lower boundary, and 
that of O 3 is fixed at 6 x 10^^ cm”^ in order to conform to measured values. 


We now turn our attention to the numerical solution of the governing equations. The method of "time 
splitting" is used to facilitate the computation; that is, the chemistry and transport portions of the 
equations are in general solved at different time steps. The technique for the "chemistry" step is dis- 
cussed first; the technique that applies to transport is described in Ref, 5., 

For altitudes of 55 km and below the chemical rate equations are solved using the following type of 
implicit technique. Equation (1) leads to the finite difference form 


1 


n. 

1 






(23) 


in which the superscript j indicates the index of the time step t. Equation (23) can be linearized 
by taking the first term in a Taylor series expansion about : 




At 


= Q. 


3n, 


(n- 


j+i 


k=i 


(24) 


The m members of the set of mass conservation equations are thus coupled and require solution by inver- 
sion of a large matrix at each grid point. Diurnal averaging of photodissociation rates and reaction 
rate coefficients is performed in a manner similar to that of the 1-D model. 

At altitudes above 55 km, the kinetic equations become sufficiently stiff that the foregoing technique 
fails. In order to solve the equations for high altitudes, we simplify the system by eliminating all 
constituents except 0(^P), O 3 , 0(^D), OH, HO 2 , H, H 2 O, and NO. To solve the kinetic equations, we 
assume 0(^P), 0(^D), and O 3 to be in photochemical equilibrium with each other, and also assume that OH, 
H 02 » and H are in photochemical equilibrium with each other. We then compute total odd hydrogen 
(HOx = OH + HO 2 + H in this case) and odd oxygen (0^ = 0(^P) + O 3 + 0(^D)) concentrations by solving the 
coupled equations for 0^ and HO^. The computed concentrations are then partitioned among member species 
by assuming photochemical equilibrium abundance ratios. 

To assess the influence of ozone on the temperature structure (and thus the chemical composition) of the 
upper stratosphere and lower mesosphere, we have included a heating and cooling code. The heating algo- 
rithm is developed from the work of Lacis and Hansen (Ref. 23). Their method, a parametric treatment 
based on accurate multiple-scattering computations, includes the effects of Earth's albedo. (Earth's 
albedo is assumed to be constant and equal to 0.34; see Ref. 24, p. 2-2). The vertical distributions of 
ozone both above and below the calculation point are accounted for in the calculation of the absorbed 
and scattered solar radiation. Because the heating rate is dependent on the zenith angle and the number 
of hours of daylight, dlurnally averaged heating rates are calculated as described in Ref. 9. The heat- 
ing rates are calculated at each latitude and altitude every 7 days so that the changing solar position 
and number of hours of daylight are taken into account. 
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The cooling rates are computed from a Newtonian cooling model (Ref. 25) such that 


R(e,p) = A(e,p) + ^ B(p)AT(e,p) . 
i 


( 25 ) 


where R is the cooling rate (in degrees/day) , AT is the model-calculated temperature minus the stan- 
dard atmosphere reference temperature, and A(0,P) and B(P) are latitude- and pressure-dependent con- 
stants. The A(0,P) were determined by requiring the radiative cooling to exactly balance the radiative 
heating for the standard reference temperature (AT = 0) and ozone profiles. The B(P) were taken from 
Ref. 25. 


In Fig. A, 2-D model predictions for the ambient 
atmosphere (autumn, lat. 40 ®N) are shown. The 
agreement with observational data is about the 
same as that obtained with the 1-D model, except 
for NO; the mesospheric concentrations appear to 
be about an order of magnitude too small between 
55 and 85 km. We attribute the low predicted 
values to the absence of a mesospheric source of 
NO and to slow upward transport of NO from the 
stratosphere; the vertical eddy diffusivity in 
the mesosphere is probably a factor of 3 too low, 
a state of affairs that is not readily correct- 
able, because of limitations in calculational 
stability. 

2.3 Simplified Plume Model 

To complement our one-dimensional (1-D) and two- 
dimensional analysis of the widespread photo- 
chemical effects caused by rocket emissions, we 
utilize a simple model of an expanding rocket 
plume to study potentially important local 
effects of the rocket plume. The model is 
applied for both rocket launch and reentry 
events, and is predicated on several approxima- 
tions. The simulated launch-reentry plume is oriented vertically and has cylindrical symmetry. The 1-D 
model is used to sample the entire plume, although the plume may be truncated in altitude extent to simu- 
late a slanted flight trajectory. The plume is given an initial width corresponding to a short period 
of expansion. The rocket injectants at each height are averaged over the initial area of the trail, and 
these concentrations, when added to the ambient concentrations, represent the starting conditions for 
the 1-D model. 

The area of the rocket trail is assumed to increase with time according to the relation (Ref. 48 )» 

A = 2K(t + to) , (26) 

where A is the area (cm^), K is the effective horizontal diffusion coefficient (cm^/sec) and time 
t = 0 at the initial expansion time of the plume (sec), to* Obviously the initial area and time can be 
related by Aq = 2Kto- The constituents under study are always taken to be uniformly distributed across 
the plume. 



CONCENTRATION, 


Fig. 4. Vertical profiles of H, OH, H 2 O, 0(^P), O 3 , 
and NO concentrations calculated with the two- 
dimensional model. The measured concentrations are 
as shown in Fig, 3. 


For a chemically active species in an expanding plume, it is easy to show that the rate of change of its 
average plume concentration, n, due to expansion only is. 


dn 

dt 


expansion 


-(n 


. 1 dA 
A dt 


(27) 


where n^ is the ambient concentration. Equation (27) presumes that chemistry occurs uniformly across 
the plume, that the plume expands by mixing in air containing an ambient abundance of each gas, and that 
mixing within the plume is instantaneous. If K is constant, equation (27), by virtue of equation (26), 
becomes , 
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dn 

dt 


n - n, 


( 28 ) 


a 


expansion ° 

This term is readily incorporated in a 1— D model as a pseudochemical process. 

The expanding cylindrical plume model just described is functionally similar to a Flckian diffusion plume 
model. In the latter case, one may consider the rate at which molecules diffuse out of an imaginary 
cylinder centered in a growing Gaussian plume. At times exceeding the initial spreading time, t , the 
two plume treatments give essentially the same solution. ° 


3. A NOCTILUCENT CLOUD MODEL AND ROCKET CONTRAILS 

With the deposition of large quantities of water vapor in the mesosphere by SPS rockets, it is likely 
that exhaust contrails will form. An important question is whether these contrails can evolve into, or 
trigger, widespread and persistent high-altitude clouds or enhance naturally occurring clouds. At high 
latitudes in the summer , the meteorological conditions at the mesopause are such that tenuous night- 
luminous (noctilucent) clouds (NLC) are a not infrequent sight. However, at low latitudes, NLC are never 
seen, except occasionally following rocket launches. 


Obviously, it would be most informative to test for the occurrence of rocket clouds from a liquid-fueled 
rocket with an appropriately designed atmospheric experiment. It is also useful, however, to make 
detailed physical simulations of rocket contrails and clouds. We have constructed a comprehensive noc- 
tilucent cloud model as part of this assessment. A summary description of the model is given in Table 3, 
and a schematic outline of the interactions between the four basic elements of the model, water vapor, 
ice crystals, meteoric dust and air ions, is given in Fig, 5. Details of the physics and photochemistry 
used in the model are given in Ref. 28 and in references therein. Reference 15 also summarizes existing 
noctilucent cloud observations and reviews previous NLC modeling work. By performing a variety of sen- 
sitivity studies with the model. Ref. 28 concludes that the model yields a reasonable simulation of typi- 
cal noctilucent cloud properties. Here, the model is used to estimate the effects of rocket water vapor 
emissions on local and global cloudiness. 

Table 3. Description of the Noctilucent Cloud Model 


Basic structure 


Basic elements 
Water vapor 


Ice crystals 


Meteoric dust 


One-dimensional model, 70-105 km in 0.5-km steps, thirty-five particle size bins from 
0.001 ym to 2.56 ym radius. Particle volume doubles between size bins. Temperature, 
density, and diffusion coefficient profiles are fixed. Solar fluxes are calculated at 
each height in 96 wavelength bins. Time-dependent or steady-state simulations. 


Fixed concentration at the lower boundary, source from methane decomposition, loss by 
photolysis. Water condenses on and evaporates from ice crystals. 

Formed by nucleation of meteoric dust, ion nucleation, and homogeneous nucleatlon, all 
involving water vapor. Microphysics includes growth, coagulation (with other ice 
crystals and dust), sedimentation, diffusion, and charging. Crystal shape effects 
(cubic, hexagonal, and amorphous ice) are also included. 

Produced from ablating meteors as 10-A smoke. Microphysics includes nucleation into 
ice crystals, coagulation, release from evaporating crystals, sedimentation, diffusion, 
and charging. 


Ionization 


Formed by galactic cosmic ray bombardment and NO photoionization. Lost by ion-electron 
and ion-ion recombination, nucleation, and collection by particles. Electrons are 
produced by photoemission. 


Figure 6 shows the predicted size distributions for ice and dust particles in an ambient noctilucent 
cloud. Also shown is the change in the ice particle size spectrum at 82 km resulting from a doubling of 
the ambient water vapor. The optical depth of the ambient (model) cloud at 550 nm is about 7 x 10"^, 
in line with values estimated from NLC observations. The optical depth of the perturbed cloud is 
1.3 X 10 . Thus, the opacity is roughly doubled when water vapor is doubled. 
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Fig. 5. Schematic outline of the noctilucent cloud model showing the major elements of the model (water 
vapor, ice crystals, meteoric dust and air ions), 'and their basic interactions. 



Fig. 6, Predicted ice crystal and meteoric dust size 
distributions in the noctilucent cloud model. Size 
distributions are given at several altitudes in the 
cloud region. The size distribution of Tozer and 
Beeson (Ref. 27) is compared to the model prediction 
(their optical measurements did not yield a unique 
size profile) . The crosshatched area indicates the 
change in the ice particle size distribution caused 
by a doubling of the water vapor concentration. 


4. CALCULATION OF NO^ PRODUCTION DURING 
HLLV REENTRY 

According to Ref. 29, the average mass of the 
HLLV at reentry is 452.6 metric tonnes. The total 
projected area of the craft is determined from 
the figures given in Ref. 29 as approximately 
2,280 m^. The initial orbit height of the craft 
is given as 450 km. Reference 29 specifies also 
that the deceleration of the HLLV during reentry 
shall be less than 1.4 g. The report does not 
specify, however, the entry flight trajectory or 
the lift and drag characteristics of the gliding 
body needed in a trajectory calculation. Presum- 
ably such aerodynamic data have not yet been gen- 
erated. In order to estimate the amounts and 
spatial distributions of NO^ produced during the 
reentry flight, therefore, the following assump- 
tions have been made: 

1. The functional relationship between the drag 
coefficient and the angle of attack of the craft 
is the same as for the Space Shuttle Orbiter, as 
given by Ref. 30. 

2. The lift coefficients of the vehicle are 
expressible as a product of the lift coefficient 
values of the Space Shuttle Orbiter given by 
Ref. 30 and an arbitrary fixed constant not 
greater than unity. Controlling lift coefficients 
in this manner would be possible, for instance, 
by varying the angle of wing flap. 

3. The angle of attack remains fixed at a high 
value initially. When the craft reaches a 
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certain Intermediate altitude (of about 70 to 80 km), the angle of attack is changed at a steady rate to 
a low value. 

4. The orbit of the HLLV is inclined 45° to the equator. 

5. The end point of the entry flight defined as the point where altitude is 50 km) is lat. 40°N and 
long. 120°W, so that the craft could land on an airfield in central California. 

6. Atmosphere rotates with the Earth without any slip. 


The entry trajectory can be varied by varying (1) the duration of the retro-rocket burn at the orbit 
(the amount of fuel required for this purpose would be very small, and so its mass is ignored, (2) the 
lift coefficient (i.e., flap angle), (3) the initial angle of attack, (4) the final angle of attack, 
(5) the altitude at which the angle of attack 


begins to decrease, and (6) the rate at which 
the angle of attack is changed. In the pres- 
ent work, these variable parameters were 
changed over a wide range in order to generate 
a family of possible entry trajectories that 
satisfies the condition of maximum allowed 
deceleration of 1.4 g. The calculation requires 
a straightforward numerical integration of the 
Newtonian equations of motion in two dimensions. 
Earth was assumed to be a perfect sphere, and 
the atmospheric data of the U.S. Standard 
Atmosphere (Ref. 31) were used for the calcu- 
lation. Figure 7 shows the profiles of such 
flight trajectories. As indicated, the trajec- 
tories are produced by varying retro-AV, that 
is, the velocity decrement achieved by firing 



the retro-rockets, between -0.105 and 

-0.110 km/s At the altitude of 120 km, these Fig. 7. Possible entry trajectories, 

two AV values produce entry velocities of 

7.588 and 7.583 km/s and entry angles of -0,896° and -1.103°, respectively. Of these, three typical 

trajectories were selected for computation of the NO^^ production rate. According to the distance 

required for the reentry flight to end, they are referred to here as "shallow,” "nominal," and "steep" 
trajectories, as shown. The Mercator projection of the entry trajectories is shown in Fig. 8. 


The rate of NO^ production along these trajectories was calculated by the method of Park (Ref. 32). In 
brief, the method assumes that the amounts of NO^ produced by a spacecraft are the same as that produced 
by a circular cone having the same overall chord length, overall surface area, and angle of attack. 

Such a method is shown in Ref, 32 to yield approximately the same amount of NO^ as a real spacecraft 
with a blunt nose, at least up to the point where such an "exact" calculation was possible. The com- 
puter code used by Park (Ref. 32) was run in the present study. The code calculates the thermodynamic 
variables simultaneously with the chemical rate equations to compute the concentrations of N and NO in 
the flow around a given cone. The computation is carried out until the flow leaves the trailing edge 
of the spacecraft. Then the assumption is made that the NO molecules are chemically frozen, but that 
all N atoms convert to NO in the wake. This assumption slightly overestimates the total NOj^ concentra- 
tion in the wake because it ignores the mutual destruction process N + NO N 2 + 0. But such error 
would be small because (1) the total concentration of N is quite small, and (2) there are far more O 2 
than NO molecules in the wake for N atoms to react with. 


The calculation yields the number of NO^^ molecules deposited along the wake per unit length of the wake 
for each entry flight. By integrating the number of the molecules along the trajectory, one obtains 
the integrated (regionally or globally) number of NO^ molecules produced. The overall efficiency of 
NOx formation is calculated from the relation 


total chemical energy of NO^ produced 


initial kinetic energy of HLLV 
total number of NO^ molecules produced x e 



( 29 ) 
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LONGITUDE 


Fig. 8. Mercator plot of the entry trajectories 
considered. 


where e is the formation energy of one NO 
molecule, e = 1.515 x 10“^^ erg, (0.95 eV) , M 
is the mass of HLLV, and is the entry 

velocity. To be consistent with the two- 
dimensional atmospheric model adopted in the 
present work, the global atmosphere was parti- 
tioned by the latitude increment of 5® and alti- 
tude increment of 2.5 km. Figure 9 shows the 
distribution of NO^ by altitude. The total 
amounts of NO^^ produced per flight, in terms of 
mass equivalent of NO, and the efficiency of con- 
version of the initial kinetic energy into NO 
are: 

Trajectory Mass of NO^, tonnes Efficiency, n 


Shallow 53.0 0.0124 
Nominal 65.5 0.0153 
Steep 71.2 0.0167 



The efficiency values obtained here are nearly 
equal to those for the Space Shuttle Orbiter 
(Ref. 32) and for large meteoroids (Ref. 33). 

Figure 10 shows the production of NO molecules 
cm as a function of altitude and latitude for 
the nominal entry case (45® angle of attack). 



Fig. 9. Vertical distributions of NO^ produced by Fig. 10. Production of NO as a function of altitude 
each HLLV reentry. and latitude for the nominal entry case (45® angle 

of attack). 
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5. A MODEL OF THE LOWER IONOSPHERE 


To assess the effect of compositional changes (l.e. , in H 2 O and NO) due to HLLV*s on the structure of 
the lower ionosphere, it is necessary to construct a model ■ that includes the appropriate ion chemistry. 
The model used, here considers a total of 24 charged species (i.e,, 0"^, ot, NO*^, N^, ot‘H 20 , H 30 “^« 0 H, 

H •H20,_H*^-(H20)2^ H‘^-^H 20 ) 3 , H+«(H 20 ) 4 , N 0 +- 02 , N 0 +-N 2 , N 0 +-H 20 , N 0 +-(H 20 > 2 , N0+- (H 2 O) 3 0" , 0^, 'Oj, , 

C 07 , CO 4 , NO 2 , N 07 , N 03 -H 20 , and electrons) in a given specified background neutral atmosphere consisting 
of N 2 , 0, O 2 , O 3 , 02 ( Ag) , NO, NO 2 , and H 2 O. The concentrations of minor species were taken from the 
1-D photochemical model (see Sec. 2.1), except that the NO concentrations between 70 and 80 km were set 
equal to 10 cm as recommended by Mitra (Ref, 34) and based on available theoretical results and 
observed ionospheric data. The basic principles used in constructing lower ionospheric models were 
described in several standard texts and monographs (e.g.. Refs. 34-37). 

Solar Lyman-a and x-rays and galactic cosmic rays constitute the sources of primary n 1, ot, o"^, and No"*" 
ions and electrons. The primary ions Initiate the ion-molecule reactions listed in Table 4; eventually, 
the positive ions are neutralized by electrons and negative ions. Negative ions are formed by electron 
attachment to neutral species (mainly O 2 and O 3 ). Attachment initiates a number of negative ion reaction 
sequences (also listed in Table 4); again neutralization is the ultimate sink for negative charge. Clus- 
tering of positive ions to neutral species such as N 2 , O 2 , CO 2 , and to H 2 O in particular, is an extremely 
important process in the D-region because it is responsible for the conversion of the primary ions ot 
and NO into hydronium ions H *(H 20 )n. Accordingly, all known important clustering processes were 
included in the model. The high air density of the D-region insures that ion-molecule reactions will 
occur rapidly. In addition, the recombination of positive and negative charge occurs in the time scale 
of -1 hr. Consequently, the charged species were assumed to be in a local steady-state condition, and 
their equilibrium number densities were obtained by solving the set of coupled chemical kinetic equations 

» 0 . (30) 

where and are the production and loss rates of the ith ion. The equations are coupled because 
the production of the ith ion generally Involves many other ions, as is evident from Table 4. Simi- 
larly, the loss of the ith ion Involves not only reaction with neutrals but also charge recombination 
processes. 


Table 4. Ion Reaction Rate Coefficients — 
Rate Coefficient = AT^ exp(-a/T)^ 


Reaction 


a 

3 

+ N 2 NO*^ + N 

3.6 (-10) 

0 

-1.0 

+ 02“^ 0^ + 0 

3.4 (- 10 ) 

0 

-0.48 

+ 0 -J- NO^ + N 

1.4 (-10) 

0 

0 

+ 02"^ 02^ + N 2 

6.4 (-7) 

0 

-1.67 

+ O 2 NO"^ + NO 

1.0 (-17) 

0 

0 

+ NO ->■ NO'*’ + O 2 

6.3 (-10) 

0 

0 

+ O 2 + M ^ ot + M 

2.0 (-30) 

0 

0 

+ 0 ot + O 3 

3.0 (- 10 ) 

0 

0 

+ H 2 O ot*H20 + 02 

2.0 (-9) 

0 

0 

H 2 O + H 2 O H 30 '''« 0 H + O 2 

9.0 (-10) 

0 

0 

H 2 O + H 2 O H"*’*H 20 + products 

3.0 (-10) 

0 

0 

r-OH + H 2 O + H'''*(H20)2 + OH 

3.0 (-9) 

0 

0 

H 2 O + H 2 O + M H"''.(H20)2 + M 

2.75(-17) 

0 

-4 

(H20)2 + M H'*'* (H 2 O) + H 2 O + M 

9.6 (11) 

17100 

-5 

(H20>2 + H 2 O + M + H^*(H20>3 + M 

1.86(-17) 

0 

-4 

(H20)a + M -»■ H'''*(H20)2 + H 2 O + M 

1.95(11) 

11000 

-5 

(H20>3 + H 2 O + M .i- H+»(H20 )i, + M 

1.94 (-17) 

0 

-4 

(H20)i, + M + H+*(H20)3 + H 2 O + M 

1.36(11) 

8360 

-5 

+ H 2 O + M -► N0'*'-H20 + M 

7.9 (-17) 

0 

-4.7 

•H 2 O + H 2 O + M + NO'^* (H 2 O) + M 

4.4 (-16) 

0 

-4.7 

• (H20>2 + M + N0'*'*H20 + H 2 O + M 

1.0 (-14) 

0 

0 

•(H20>2 + H 2 O + M + N0'^.(H20)3 + M 

4.4 (-16) 

0 

-4.7 

• (H 2 O) 3 + M + NO'^* (H20>2 + H 2 O + M 

1.9 (-12) 

0 

0 
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Table 4 . Concluded 



. Reaction 


a 

3 

NO 

(H2O) 3 + H2O + (H2O) 3 + HNO2 

7.0 

(-11) 

0 

0 

NO 

^ + CO2 + M -> N0'''-C02 + M 

1 . 89 (- 22 ) 

0 

-3 

NO 

"^•C 02 + M “>■ no"*" + CO 2 + M 

3.1 

( 4 ) 

4590 

-4 

NO 

^•C02 + H2O ^ N0“''-H20 + CO2 

1.0 

(- 9 ) 

0 

0 

NO 

^ + N2 + M -> N 0 '^*N 2 + M 

2.0 

(- 31 ) 

0 

- 4.4 

NO 

+ M + NO"'' + N2 + M 

1.5 

(6) 

2450 

- 5.4 

NO 

^•N2 + CO2 + N0‘''»C02 + N2 

1.0 

(- 9 ) 

0 

0 

NO 

^•N2 + H2O N0‘^*H20 + N2 

1.0 

(- 9 ) 

0 

0 

O3 

+ e -> 0” + O2 

1.0 

(-12) 

0 

0 

0“ 

+ O3 ->■ O3 + 0 

7.0 

(-10) 

0 

0 

0“ 

+ 02+02*^03+02 

1.0 

(- 30 ) 

0 

0 

02 

+ e + N2 ^ 01 + N2 

1.0 

(- 31 ) 

0 

0 

02 

+ e + O2 O2 + O2 

1.8 

(- 30 ) 

0 

0 

01 

+ 0 *> O3 + e 

3.3 

(-10) 

0 

0 

01 

+ O2 O2 + O2 + e 

2.0 

(-10) 

0 

0 

01 

+ O3 -► O3 + O2 

3.0 

(-10) 

0 

0 

01 

+O2 +M ^04 +M 

1.0 

(- 30 ) 

0 

0 

02 

+ CO2 + 02 “^ CO4 + O2 

2.0 

(- 29 ) 

0 

0 

01 

+ O O2 + O2 

1.4 

(-10) 

0 

0 

01 

+ C02 C03 + 02 

4.0 

(-10) 

0 

0 

CO' 

3 + 0 O2 + CO2 

8.0 

(-11) 

0 

0 

CO' 

3 no -y NOI + CO2 

9.0 

(-12) 

0 

0 

CO' 

3 + NO2 ~y no's 0O2 

8.0 

(-11) 

0 

0 

NO' 

2 + O3 N07 + O2 

1.8 

(-11) 

0 

0 

04 

+ no NOi + O2 

2.5 

(-10) 

0 

0 

04 

+ CO2 ~y CO4 + O2 

4.3 

(-10) 

0 

0 

co; 

J + NO -V NO3 + CO2 

4.8 

(-11) 

0 

0 

CO4 + 0 col + O2 

1.5 

(-10) 

0 

0 

NO' 

3 + H2O + O2 -y N0T-H20 + O2 

7.5 

(- 29 ) 

0 

0 

n| 

+ e products 

2.9 

(- 7 ) 

0 

0 

02 

+ e products 

2.0 

(- 7 ) 

0 

-1 

NO 

^ + e products 

4.0 

(- 7 ) 

0 

-1 

XY' 

^ + e ->• products^ 

1.0 

(- 5 ) 

0 

0 


+ Y” -»■ products^ 

1.0 

(- 7 ) 

0 

0 

w“ 

+ hv -»■ products + e 

0.2 


0 

0 


^T denotes atmospheric temperature, 
^In units of sec"^ for unimolecular 

reactions, 

cm^ sec 

-1 


for bimolecular reactions; and cm® sec"^ for thermolecular 


reactions. 

All other dissociative recombination processes. 
^All ion-ion recombination processes. 


A fast stiff differential equation solver, due to Gear (Ref. 38 ) was used to obtain the solutions of 
equations ( 30 ) . The production and loss mechanisms of each species were determined automatically by the 
computer code by scanning the library of chemical reactions used in the model. Details of this computa- 
tional scheme are described by Prasad and Huntress (Ref. 39 ) in the context of their model of interstel- 
lar chemistry, in which 137 species and more than 1,400 chemical reactions are involved. 


Table 5 presents the values of the electron number density n^, and the ratio of the concentration of 
negative ions to the concentration of electrons, X, at various altitudes, with ambient abundances of 
water vapor and nitric oxide computed with the photochemical models. The data in Table 5 simulate 
globally averaged daytime equilibrium conditions. Our quiescent values of n^ and X are consistent with 
observations and with the results from other models as described in the previously referenced texts and 


14 



monographs. This demonstrates the adequacy 
of the model for the present assessment where 
only possible changes in the ionosphere due to 
alterations in the neutral atmosphere caused 
by HLLV operations are sought. 


Table 5. Simulated Unperturbed Daytime Globally Aver- 
aged Ionosphere 


6 . 


THE ACCUMULATION AND DISPERSION OF 
ROCKET EXHAUST CONTAMINANTS 


Altitude, 

km 

[H 20 ], 
cm ® 

[NO], 

cm”^ 

[e] 

cm” 

9 

3 

tx"] 

cm" 

a 

3 


100 

1.3 (7)° 

1.0 

(8) 

4.5 

(4) 

1.6 

(0) 

3.6(-5) 

90 

1.3 (8) 

4.8 

(7) 

7.5 

(3) 

4.7 

(0) 

6. 3 (-4) 

85 

4.1 (8) 

1.5 

(7) 

4.1 

(3) 

3.6 

(0) 

8.6(-4) 

80 

1.3 (9) 

1.0 

(7) 

8.0 

(2) 

1.0 

(1) 

1.3(-2) 

75 

3.0 (9) 

1.0 

(7) 

6.1 

(2) 

7.1 

(1) 

1.2(-1) 

70 

7.3 (9) 

1.0 

(7) 

1.7 

(2) 

9.8 

(1) 

5.8(-l) 

65 

1.5(10) 

9.6 

(6) 

2.9 

(1) 

1.2 

(2) 

4.2 (0) 

60 

3.0(10) 

6.5 

(7) 

7.6 

(0) 

2.0 

(2) 

2.6 (1) 


^[X ] E concentration of all negative ions. 


h = [X ]/[e]. 

(7) = 1.3 X 10’ 


The models described above have been used to 
simulate the accumulation and dispersion of 
the substances deposited during HLLV launch 
(water vapor, nitric oxide, carbon dioxide, 
carbon monoxide, sulfur dioxide) and created 
during reentry (nitric oxide) . In order to 
simulate the long-term effects, the 1-D and 
2-D photochemical models were run for 
10 years of continuous HLLV operations at 
the rate of 4000 launches per year with the 

injectants averaged over different horizontal extents (Northern Hemisphere, the globe, etc.). The 1-*D 
model simulations appeared to have reached a steady state. The 2-D model simulations had closely 
approached a steady state in the mesosphere, but had not quite attained steady state in the stratosphere. 
The 1-D model directly simulated water vapor deposition at all altitudes up to 120 km. Water vapor 
deposition at higher altitudes was simulated by inserting it into the model at 120 km altitude from where 
it could be transported downward or converted photochemically to H 2 and flow upward as well as downward. 
The H 2 could escape as discussed below. In the 2-D model, which did not contain H 2 , all water deposited 
above 90 km altitude was assumed to be injected at 90 km; the upper boundary condition was constructed 
such that water vapor could move across the boundary seeking equilibrium of high-altitude water vapor 
(l.e., above 90 km) with the water vapor below. The short-term effects are treated with the plume (dis- 
persion) model discussed in Sec. 2. In the present section we discuss the widespread accumulation of 
injectants under quasi-steady-state conditions and the dispersion of the water vapor contrails and 
reentry nitric oxide plumes for single events. 


6. 1 Water Vapor 

The predicted increases in the water vapor abun- 
dance obtained with the 1-D and 2-D photochemi- 
cal models for 10 years of HLLV operations are 
shown in Fig. 11. It is apparent from the 
results of the 2-D model that no significant 
"corridor effect" (enhanced accumulation in 
the vicinity of the launch latitude) is likely 
at altitudes below 70 km. At an altitude of 
85 km (the highest altitude at which the 2-D 
model yields reliable results) the computed cor- 
rodor enhancement is only -15%, which suggests 
that corridor effects are probably negligible 
at middle atmospheric levels. Thus the changes 
in air composition caused by HLLV water vapor 
deposition can be adequately simulated with 1-D 
photochemical models. The globally averaged 
increase in H 2 O obtained with the 1-D model is 
shown by the broken curve in Fig. 11. In the 
stratosphere (i.e. , at altitudes below 50 km) 
the two models predict nearly the same increase 
in water vapor. However, between 50 and 75 km 
(lower mesosphere) , the increase obtained with 
the 2-D model is about double that obtained 
with the 1-D model. The difference is a conse- 
quence of the difference in vertical transport 
rates and the fact that the 2-D model does not 
simulate conversion of water vapor to H 2 ; 
the vertical "eddy diffusion" coefficient in 


PERCENT INCREASE OF H 2 O AT 30°N 
5 10 15 



Fig. 11. Computed absolute increase in the water 
vapor concentration after 10 yr of launches of HLLV 
at a rate of 400 yr”^. The solid lines represent the 
results from the 2-D model and the broken line the 
globally averaged computed increase from the 1-D 
model. 
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the 2-D model is somewhat smaller than that in the 1-D model. On the other hand, the actual value is 
uncertain and our results give an Indication of the effects of this parameter on the H 2 O prediction. 

Many meteorologists believe that the influx of water vapor to the stratosphere is controlled by the low- 
temperature cold trap" of the tropical tropopause (Refs. 14, 40, 41). The mechanisms that remove water 
vapor from the stratosphere are still unknown, but are probably dominated by subsidence at mid-latitudes 
and high latitudes. We have made such an assumption in our model calculations. Another loss mechanism 
for water vapor is photochemical conversion to relatively inert molecular hydrogen in the mesosphere, 
mainly through the reaction sequence. 


H 2 O + hv OH + H 

(31) 

OH + O 3 ->■ HO 2 + O 2 

(32) 

H + O 2 + M HO 2 + M 

(33) 

H + HO 2 H 2 + O 2 

(34) 

Some of the H 2 diffuses to the troposphere, and some is dissociated at high altitude, contributing to 
the hydrogen escape flux to space (for discussion of the escape mechanisms, see Ref. 42). According 
to recent measurements the hydrogen escape rate is -10® H atoms cm“^ sec“^ (Ref. 43). An escape flux of 
this magnitude was Imposed as an upper boundary condition on the hydrogen in our 1-D model. 

Molecular hydrogen may also be converted back into water vapor by reactions such 

as 

0(^D) + H 2 OH + H 

(35) 

OH + HO 2 ^ H 2 O + O 2 

(36) 


However, in the mesosphere, H 2 is produced at the expense of H 2 O, and diffuses into the stratosphere and 
eventually the troposphere. 

The effect of the deposition of water vapor and hydrogen at high altitudes (up to 120 km) is a substan- 
tial thermospheric buildup of molecular and atomic hydrogen. The H 2 dominates below -120 km and the H 
above 120 km (e.g. , see Refs. 43, 44). The computed increase in total hydrogen above 100 km (globally 
averaged) is one of about a factor of 2 (i.e., a doubling of thermospheric H and H 2 ). The implication of 
such a large increase in thermospheric hydrogen worldwide are beyond the scope of the present study, but 
must be boldly underscored. It has been suggested to the authors by B, M. McCormac^ that the population 
of ring current protons in the magnetosphere could be substantially altered by charge exchange. We men- 
tion here that because combustion is not stoichiometric, hydrogen is deposited along with water vapor in 
a volumetric ratio of 0.3. The 1-D model takes account of H 2 photochemistry, but the 2-D model does not. 

Because the short-term local dispersion of water vapor rocket emissions is associated principally with 
the production of "contrails" (small ice particles), we defer discussion of this problem to Sec. 8 . 

6.2 Nitric Oxide 

Nitric oxide deposited in the stratosphere during the launch phase of the HLLV leads to an Increase in 
the global NO abundance of less than 1% over ambient concentrations for 400 launches per year. 

Mesospheric nitric oxide is not normally produced there. Rather, it is produced by photochemical pro- 
cesses in the stratosphere and the thermosphere and then transported upward or downward, respectively. 

The stratospheric source results from the reaction of 0(^D) with nitrous oxide, and thermospheric nitric 
oxide is produced through the dissociation of N 2 , mainly via ionic processes. 

However, NO created in the mesosphere and lower thermosphere during HLLV reentry can lead to substantial 
increases in its upper atmospheric abundance. The relative magnitude of the Increase depends on the 
background NO concentrations which are quite uncertain. Figure 12 shows the predicted absolute Increases 
in NO concentrations caused by 10 years of HLLV launches at the rate of 400 per year and reentry on the 
nominal trajectory. These represent less than 20% increases in NO at all heights. 


^Lockheed Palo Alto Research Laboratory, 1980. 
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The "corridor" effect for NO is much more pro- 
nounced than was the case for water vapor. 

Nitric oxide has a rather short lifetime in the 
mesosphere, about 4 days at an altitude of 
70 km and even shorter at higher altitudes. This 
is a consequence of rapid photolysis by solar 
ultraviolet radiation (refs. 7, 45), 

NO + hv N + 0 , (37) 

followed by the odd-nitrogen destruction reac- 
tion, 

N + NO N 2 + 0 . (38) 

Of course, reaction (37) must compete with the 
NO recycling reaction 

N + O 2 NO + 0 , (39) 

which inhibits NO loss; the rates of reac- 
tions (38) and (39) are comparable at the 70-km 
altitude. 

The NO lifetime can therefore be expressed as 

^ ^ k:^B[N0] + k:39[02] /40) 

^NO " 2J37k3g[N0] 

In equation (40) , k denotes a rate coefficient and J a photolysis rate; values of the reaction coeffi- 
cients were taken from Ref. 6 and photolysis rates from Ref. 7. Since the concentration of O 2 decreases 
rapidly with ascending altitude, the lifetime also decreases, reaching an asymptotic value of ~2 days at 
90 km. One can account for the relatively strong enhancement of NO in the region of rocket reentry by 
noting that the NO is likely to be destroyed by the sequence of reactions (37) and (38) before it can be 
transported raeriodlonally over a significant distance. For the transport parameter values used in the 
2-D model, the mean distance through which material is transported meridionally in 4 days is only 
"400 km. 

The mesospheric source of nitric oxide includes photochemical reactions which dissociate N 2 , both EUV 
and X-rays (Ref. 46) as well as solar protons, meteorites (Ref. 33), etc. The 1-D model includes such 
sources of NO as well as a flux from above 120 km of •'10^ molecules cm^ sec 

The NO predictions from the 1-D and 2-D models are in rough qualitative agreement, but differ in details. 
Because it is not possible to simulate the reentry trajectory in the 1-D model, substantial differences 
between the results from the two models are expected. Nitric oxide also cools the thermosphere in the 
5.3-pm emission band (Ref. 47), but the amount of NO produced above 120 km by HLLV reentry is very small, 
and this effect is expected to be negligible. 

We have also simulated the dispersion of an NO reentry trail with the aid of the plume model discussed 
in Sec. 2. Figure 13 illustrates the dispersal of NO at points along the reentry flightpath where NO is 
generated in large quantities (see Figs. 7,9). The horizontal eddy diffusivity employed for the calcu- 
lations was 10® cm^ sec"^; actually, K is scale-dependent and hence is time-dependent_(Ref , 48), its 
effective value increasing with time. For periods of less than 1 day, K ~ 10® cm^ sec ^ is a reason- 
able value, but it increases rapidly thereafter, with the result that the trail undergoes rapid dispersal 
to concentrations approaching ambient values. The effect of such large NO enhancements on ionospheric 
structure will be discussed later. 

6.3 Carbon Dioxide and Carbon Monoxide 



The first stage of the HLLV burns methane fuel and carbon dioxide is a principal product of combustion. 
Since the methane molecule contains carbon and hydrogen atoms in the ratio of 1 to 4 , assuming stoichio' 
metric combustion one expects the volume rate of CO 2 emission to be one half that of water vapor at 
altitudes below 56 km. Actually, according to Ref. 1, the combustion is not stoichiometric; the CO 2 
volume emission rate is about 0.4 that of water vapor and that for CO is about 0.3 that for water 
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NO CONCENTRATION, cm-3 


Fig. 13. Nitric oxide concentrations in a dispers- 
ing trail formed during the reentry of an HLLV, The 
plume model described in the text was used. The 
horizontal diffusion coefficient is 10® cra^ sec“^. 
Three sections of the plume along the reentry tra- 
jectory are shown. The ambient 1-D model nitric 


(Ref. 1). Although CO2 deposition ceases at an 
^ll^itude of 56 km, it is effectively mixed upward 
from that height because its only sink, diffusion 
to the troposphere, is a slow process. As a 
result, the absolute increase in the CO2 content 
should roughly approximate one half that for water 
vapor. Figure 14 shows the computed change in the 
mixing ratio of CO2 for several altitudes. The 
increase in the mixing ratio of CO2 at an altitude 
of 50 km therefore ranges from about 3 x io“® ppmv 
to about 7 X 10 ^ ppmv, which should be compared 
with the mixing ratio in the ambient atmosphere 
of "300 ppmv. The fractional change is only about 
1 to 2 parts in 10®, which seems to be much too 
small to cause a significant increase in the 
stratospheric-mesospheric cooling rate. If all of 
the cooling in this region were due to emission of 
infrared radiation by CO2 » a very rough estimate 
of the resulting temperature change is "lO”® K or 
less. The heating-cooling models described in 
Sec. 2.1 were used to confirm this estimate. 

Carbon monoxide emitted by HLLV’ s cannot accumu- 
late like CO2 because CO is rapidly oxidized to 
C02« Afterburning in the rocket plume should 
convert most of the emitted CO into CO2 . ^ Other- 
wise, the CO will mix into the atmosphere and be 
oxidized principally by reacting with OH, 

CO + OH CO2 + H , (41) 


oxide profile is also given. with a rate coefficient ~1 x 10“^® cm® sec“^, 

corresponding to an upper atmospheric lifetime of 
^CO < 3 months. Reactions of CO with O2, O3 , 

HO2 , etc., are extremely slow and can be Ignored. 
Inasmuch as stratospheric OH concentrations are 
-1 X 10^ molecules cm“® (Ref. 13), the lifetime 
of injected CO could be as short as 2 weeks, 
assuming OH levels are not greatly perturbed. CO 
concentration Increases would have to exceed 
lO®"® molecules cm“® to affect OH concentrations 
(background amounts of CO in the stratosphere 
vary from -2 x 10^® molecules cm”® near 20 km to 
'*1 X 10® molecules cm”® near 40 km based on 
recent model calculations). Thus, if significant 
amounts of CO remained in the rocket plume, the 
_-Q o OH abundance near the plume could be reduced by 

"I® 10" 10" reaction (41). However, OH would also be influ- 

CHAIMGE IN MIXING FRACTION OF CO 2 enced by the large quantity of water contained in 

the plume, which might compensate the CO effect. 

Fig. 14. Computed increase in the mixing ratio of 

carbon dioxide after 10 years of HLLV launches at In the near-wake of the rocket, and in the launch 

a rate of 400 per year. corridor, CO buildups might be larger. High con- 

centrations of emitted CO could reduce OH levels, 

lengthening the residence time of CO in the exhaust plume. However, the plume would be rapidly diluted 
with ambient air as it disperses, allowing the OH reaction cycle to begin depleting CO quite efficiently, 
perhaps at a rate approaching (or, locally, exceeding) the global oxidation rate(s) estimated above. In 
addition, because actual CO emissions are probably much lower than the upper-limit values just consid- 
ered, CO perturbations by SPS rockets will very likely be negligible. Calculation with the 1-D model 
for the globally averaged case indicates that the CO abundance would be increased by less than 10%, which 
would have a negligible effect on the distributions of other constituents. 



®Based on calculations supplied by K. L. Brubaker, Argonne National Laboratory. 
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6.4 Sulfur Dioxide 


It has been estimated** that the first-stage fuel of the HLLV could contain 0.05% sulfur by mass. The 
yearly deposition of sulfur dioxide in the stratosphere is therefore estimated to be not more than 
"200 tonnes yr~^ , as compared with an estimated ~2 x 10^ tonnes yr ^ of SO2 produced in the stratosphere 
from COS of tropospheric origin (Ref. 49). Based on the work reported in Ref. 49, we estimate that the 
resulting fractional increase in the mass concentration of the stratospheric aerosol layer would be 
negligible and not affect Earth’s albedo and mean surface temperature. 


7. PERTURBATIONS TO ODD-OXYGEN AND ODD-HYDROGEN 


The computed accumulations of excess water vapor and nitric oxide lead to predicted changes in the dis- 
tribution of odd-oxygen and odd-hydrogen (HOx) species. The odd-hydrogen concentration is enhanced 
mainly through the reactions 

0(^D) + H2O ^ 20H (^2) 

H2O + hv ^ OH + H (31) 


where water vapor photolysis occurs at high altitudes (above 70 km) . 


The odd-hydrogen generated is distributed among 
several species including H, OH, H02» and, at low 
altitudes, H2O2 and HNO3, Our computations show 
that the increases in HOx abundances are signifi- 
cant only in the mesosphere and thermosphere; in 
the stratosphere, fractional OH and HO2 increase 
only about 0.5% at 50 km, and decrease with 
descending altitude. Of course such strato- 
spheric perturbations will influence the ozone 
column density somewhat, as is discussed below. 
Figure 15 shows computed absolute increases in 
HOx species concentrations. The calculations 
were carried out with the 1-D model for globally 
averaged conditions, which are approximate 
according to the 2-D results in Fig. 11, indi- 
cating no strong corridor effect. Although the 
results for OH and HO2 from the 1-D and 2-D 
model agree quite well up to 60 km, they diverge 
at higher altitudes. The predictions made with 
the 1-D model are probably more reliable because 
the chemical scheme used was more complete than 
the scheme used in the 2-D model, H2 was included 
water vapor deposition above 90 km more realistic, 
dictions are not expected to be reliable in any ca 



Fig. 15. Predicted absolute changes in H, OH, and 
HO2 concentrations due to water vapor emissions by 
HLLV’s. The results correspond to 10 yr of emissions 
for 400 launches yr“^. 

i the 1-D model, and the technique for simulating 
At altitudes above 80 to 85 km, the 2-D model pre- 
b because of proximity to the upper boundary. 


In the region near the mesopause, excess injected atomic hydrogen may influence the thermal balance. 

The chemiluminescent reaction with ozone, 

H + O3 OH* + O2 , (^3) 


would enhance hydroxyl (Meinel) vibrational band emission and cool the atmosphere. More importantly, 
hydrogen catalyzes odd-oxygen recombination through reaction (43) followed by 

OH + 0 H + O2 . (^^) 

The net reaction cycle releases as heat the chemical energy stored in atomic oxygen, providing an addi- 
tional local heating source. We have not been able to quantitatively assess the Importance of this heat 
source because of complications involving turbulent heat transport. 


**By K. L. Brubaker, Argonne National Laboratory, 1980. 
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In the mesosphere and upper stratosphere, odd-oxygen species, O 3 , 0(^P), and 0(^D), are strongly influ- 
enced by the catalytic reaction sequences (43) and (44), in addition to the Chapman oxygen reactions. 
Other important catalytic cycles are 


H + O 2 + M HO 2 + M (33) 

OH + 0 H + O2 ( 44 ) 

HO2 + 0 OH + O2 ( 45 ) 

net: 0 + 0 O 2 ; 

and 

OH + O3 ->• HO2 + O2 ( 32 ) 

HO2 + 0 OH + O2 ( 45 ) 

net: 0 + O3 2O2 . 


The stratosphere is more complex than the mesosphere, but several key odd-oxygen reaction cycles have 
been identified: 

OH + O 3 “>■ HO 2 "f" O 2 (32) 

HO 2 + O 3 OH + 2 O 2 (46) 


net: 03 +03 ->- 3 O 2 ; 

and 

NO + O 3 NO 2 + O 2 (47) 

NO 2 + 0 NO + O 2 (48) 


net: 0 + O 3 2 O 2 . 

Interestingly, nitric oxide deposited in the lower stratosphere may lead to ozone increases (Ref. 50) 
through the decomposition of peroxy compounds, e.g.. 


NO + HO 2 

^ NO 2 + OH 

(49) 

NO + CH 3 O 2 

NO 2 + CH 3 O 

(50) 

NO 2 + hv 

NO + 0 

(51) 


The critical reaction (49), also Interferes with the odd-oxygen destructive cycles (32) and (45), and 
(32) and (46). 

We have computed ozone column changes due to water vapor and nitric oxide Injection by HLLV launch and 
reentry operations. The computed globally averaged ozone reductions due to water vapor deposition 
(400 launches yr“^) are about 0.01%, with no apparent corridor effect. This is less than one tenth the 
ozone reduction computed for the Space Shuttle.^ Similar calculations have been made for nitric oxide 
deposition during atmospheric reentry. On the other hand, the predicted increase in ozone obtained with 
the 1-D model is only '•0.01%. We attribute the difference between the two results mainly to differences 
in the simulation of transport; the 2-D model simulates horizontal transport, whereas the 1-D model does 
not. Of lesser importance are differences in mean sun angle (the 1-D model uses a mean solar zenith 
angle which is time invariant, while in the 2-D model, the solar zenith angle is latitude- and time- 
independent) and to differences in lower boundary conditions. In the 2-D model, NO is transported 


^R. C. Whitten, unpublished computations. 
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downward into the stratosphere and then stored below altitudes where it is photodissociated . Because of 
reaction (49) an ozone increase is predicted, about 0.02% globally averaged, using the 2-D model. The 
computed latitude dependence of the column increase is shown in Fig. 16a for two seasons. It should be 


noted that such increases are highly uncertain, 
and it may be just as likely that small ozone 
decreases would occur for NO injection in the 
lower stratosphere (Ref. 51). Figure 16b shows 
predicted changes in ozone concentration as a 
function of altitude. 

Figure 17 shows the computed changes in the 
ozone and atomic oxygen concentrations above 
40 km. There is reasonably good agreement 
between the globally averaged 1-D model results 
and the 2-D model results. The predicted changes 
in mesospheric odd oxygen species do not seem to 
be especially significant, except perhaps with 
respect to chemical heating near the mesopause; 
reductions in 0(^P) at the mesopause may lead to 
slightly lower mesospheric temperatures, but the 
magnitude of the cooling, although difficult to 
estimate, is probably small (recall that the 
0 concentration is smaller because the H concen- 
tration is larger, so that the overall rate of 
the exothermic reaction cycle (43) and (44) is 
not altered as significantly). 



(a) Predicted relative increase in the total column 
ozone density as a function of latitude and season 
corresponding to NO^^ production due to reentry of 
400 HLLV*s yr“^. 


8. ROCKET CONTRAILS AND NOCTILUCENT CLOUDS 

Meinel et al. (Ref. 52) reported seeing meso- 
spheric clouds to the west of Tucson, Arizona, 
on several evenings in 1963. The sightings were 
unusual in that natural NLC*s are never observed 
at such low latitudes. Meinel et al. correlated 
the events with rocket launches in southern Cali- 
fornia, and attributed the clouds to water vapor 
and particulates deposited near the mesopause by 
rocket engines. They likened the clouds to NLC. 

Benech and Dessens (Ref. 53) also observed arti- 
ficial noctilucent clouds (ANLC) on two evenings 
when test rockets were being accelerated through 
the upper mesosphere. “On the second occasion, 
careful observations were made of the evolution 
of the cloud , which was visible against the 
darkening twilight sky for more than 30 min. 
Initially, a bright exhaust trail appeared, but 
this quickly dissipated except for a short seg- 
ment (~10 km in length) near the mesopause height. 
The segment elongated over a 50-km horizontal 
distance and eventually became distorted. 

Benech and Dessens state that the rocket deposited 
260 g H 2 O, 520 g HCl, and 1,600 g AI 2 O 3 along 
the 10-km segment of the trail that was persis- 
tent. It is interesting to note that 1 kg of 
H 2 O can be dispersed into 0.1 pm ice spheres to 
produce a 30 km^ cloud with a vertical optical 
thickness of about 10“**; that is, a small nocti- 
lucent cloud. The 1,600 g of AI 2 O 3 alone would 
not be sufficient to form a visible cloud, 
because the ejected particles are typically much 
smaller than 0.1 pm, and are optically transparent 



(b) Predicted changes in ozone concentration obtained 
with the 1-D and 2-D models. 

Fig. 16. Ozone increases. 
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(Ref. 54). On the other hand, AI2O3 particles are 
excellent ice sublimation nuclei. The combination 
of H 2 O and HCl could form a condensate of lower 
vapor pressure than pure water. The supersatura- 
tion required to nucleate 0.01 ym AI 2 O 3 particles 
would also be lower than that required to nucleate 
the 0 . 001 -ym meteoric dust normally found in the 
mesosphere, Benech and Dessens suggest in addi- 
tion, that after sunset, the AI2O3 particles may 
be cooler than the ambient air due to efficient 
radiative cooling; this would give AI 2 O 3 a 
further advantage as an ice nucleating agent. 

Past observations of mesospheric clouds associated 
with rocket launches, therefore, indicate that 
tenuous short-lived clouds are generated by solid- 
fueled rockets, which emit copious quantities of 
sublimation nuclei. To our knowledge, mesospheric 
clouds have never been observed following any of 
the large Saturn V liquid fuel (H 2 -O 2 ) rocket 
launches from Cape Kennedy, Florida. The major 
exhaust constituent of those rockets is H 2 O. 

). Admittedly, launches from Cape Kennedy could pro- 
duce clouds far to the east of Florida where the 
prevailing winds could dissipate them before the first twilight. Yet, no reports of unusual night- 
luminous clouds, which might be expected with such large water injections, have been uncovered. Accord- 
ingly, there is little evidence linking large injections of (pure) water into the mesosphere with exten- 
sive persistent, optically-thick clouds. 



Fig. 17. Predicted changes in ozone and 0(^P) due 
to water vapor emission by HLLV’s (400 launches yr 


Based on the projected SPS water vapor increases reported in Sec. 5, we expect less than a factor of 
2 increase in water vapor between 80 and 90 km over the Northern Hemisphere due to SPS rocketry. Hence, 
the likely increase in the (average) noctilucent cloud opacity is <10“^. We estimate a corresponding 
surface temperature change <0.01 K (Refs, 55, 56), which is negligible. 

The southern extent in latitude of noctilucent clouds, likewise, should not be greatly affected by SPS 
rockets. Normally, these clouds occur sporadically at latitudes as low as 45'’. Theon et al. (Ref. 57) 
noted that NLC do not automatically occur whenever the mesosphere is extremely cold. They proposed that 
another factor, probably either an enhanced water vapor concentration or a suitable condensation nucleus 
is needed to form NLC. Thus, the change in water vapor induced by SPS rockets might Increase slightly 
the frequency of NLC below 70®, without a measurable climatic effect. Above 70® latitude, Donahue et al. 
(Ref, 58) detected (from the OGO -6 satellite) a persistent circumpolar cloud layer at 80 km in summer. 
Inasmuch as models of atmospheric circulation indicate substantial upward convection of air at the summer 
poles (Ref. 59), the likely source of water for the circumpolar clouds is the polar upper stratosphere. 

In Sec. 6 it was found that SPS rockets will not noticeably disturb the stratospheric water vapor reser- 
voir. Thus, rocket H 2 O emissions should have an even smaller effect on circumpolar clouds than on low- 
latitude NLCs. 


Using the rocket plume "expansion" model outlined in Sec. 2.3 and the noctilucent cloud raicrophysics 
model described in Table 2, a rough simulation of an HLLV launch plume was made. We Ignore the complex 
aerodynamical effects in the early development of the rocket trail, and assume that after a short relax- 
ation time (-100 sec) the trail has reached temperature and pressure equilibrium with the ambient atmo- 
sphere and has a cross-sectional area of -1 km^ At this time, the exhaust water vapor is uniformly 
distributed across the plume, which is centered on the vertical axis of the 1-D NLC model. Of course, 
the initial concentration of water vapor in the plume is not allowed to exceed the local air density. * 
Normally, contrails are seen just behind rockets, so that our initial spreading time of -100 sec is 
somewhat artificial; on the other hand, the initial adiabatic expansion of the plume is very rapid and 
cooling and nucleation can occur earlier. 

In the plume model, a burst of homogeneous water nucleation occurs, and, because of the high supersatura- 
tion, these particles grow rapidly to sizes of -1 pm (radius). As a result, the water vapor in the plume 
is consumed, the supersaturation over the (presumed) frozen particles falls to - 0 , and further nucleation 
and growth ceases. It is important to note that the vapor pressure over ice is significantly lower than 
the vapor pressure over liquid water. Homogeneous nucleation produces supercooled water droplets, and. 
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because the supercooled droplets may not freeze immediately, they may remain liquid throughout the growth 
phase. As the droplets froze, the contrail might take on the appearance of a cirrus cloud. 

Subsequent to the termination of particle growth, the dominant contrail process is the expansion, which 
rapidly dilutes the particle concentration. Interestingly, the water vapor concentration is not diluted, 
because the water vapor remains in equilibrium with the particles, which evaporate and maintain a con- 
stant concentration of the vapor phase. Coagulation and sedimentation are important secondary processes. 
For example, the sedimentation rate of 1-ym ice spheres at 80 km is about 60 km day"^. We have neglected 
temperature feedback effects in the plume (infrared and ultraviolet radiation, photochemistry, and latent 
heat); proper treatment of these effects is beyond the scope of this assignment. 

As the rocket trail continues to expand, the ratio of gaseous to condensed water increases steadily at 
each height (for a fixed ambient temperature) . Eventually a point is reached where the particles begin 
to evaporate rapidly to maintain the vapor pressure of the water. At a given altitude this occurs when 
the plume has expanded to such a size that the total water vapor injected, if uniformly distributed over 
the plume area, would just produce ice saturation. Hence, after a short time, the contrail is confined 
to the coldest region of the mesosphere (assuming that the rate of horizontal expansion is roughly con- 
stant with height). The contrail also disappears abruptly when the water vapor in the coldest region 
becomes sufficiently diluted (allowing that the ambient water vapor concentration is moderately under- 
saturated). No persistent clouds are predicted because pure ice crystals cannot form and grow in the 
absence of sufficient water vapor saturation. Meteoric dust is present in the mesosphere (Ref. 60), but 
is probably too small (<0.005 yra) to act as sublimation nuclei unless the supersaturation exceeds '-lO 
(Ref. 28). 

Figure 18 summarizes the HLLV plume calculation. 

The initial optical depth along the axis of the 
contrail exceeds 100. The opacity falls off 
rapidly as the plume expands. To be clearly 
visible in daylight, the opacity through the 
width of the contrail must exceed 0.1. Note 
that, while the opacity along the axis of the 
trail decreases rapidly with time (approximately 
as A”^, where A is the cross-sectional area), 
the opacity across the trail decreases more 
slowly (roughly as A”^'^), Our calculations 
indicate that the trail could be visible for 
12 hr and cover an area of several hundred square 
kilometers. The actual horizontal extent of the 
contrail could be much larger because of the 
slant of the launch trajectory and the stretching 
and twisting of the trail by wind shears. 


After several hours, the opacity of the trail 
begins to decrease at a faster rate because the 
ice particles at some altitudes are spontaneously 
evaporating. 

A subvisible cloud could persist for about 1 day 
before the water vapor saturation point was 
reached at all heights. This remnant contrail 
might produce a spectacular noctilucent cloud 
display at first dust or dawn after an HLLV 
launch. 

In summary, we find that the effect of SPS rocket water vapor emissions on global cloud coverage is 
negligible from the point of view of climatic change. Preliminary simulations of rocket plumes indicate 
that contrails will form behind HLLVs and that they might persist for day. However, the contrails 
would cover only a small global area and have little climatic significance. Except for the possibility 
of impressive local noctilucent cloud displays in the launch corridor, we have uncovered no evidence for 
sustained, widespread clouds caused by liquid propellant rockets. We stress, however, that our conclu- 
sions are based on limited scientific data and approximate computer models and should be considered pre- 
liminary and tentative. 



Fig. 18. Simulated HLLV launch plume characteris- 
tics. Shown are the (vertical) optical depth along 
the axis of the plume at 550 nm, and the cross- 
sectional area of the plume as a function of time 
after passage of the rocket. The average size of 
the ice particles in the plume and the height extent 
of the contrail are indicated in the figure. 
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9. PERTURBATIONS OF THE LOWER IONOSPHERE 


9 . 1 Long-Term Global Effects 

Increased water vapor in the ionospheric D-region due to HLLV launches is expected to increase the rela- 
tive abundance of positive ions clustered to water molecules and to decrease the electron concentration, 
ne, because the dissociative recombination of cluster ions is very fast. Our results, however, show no 
significant changes either in the electron density profile or in the total number densities of cluster 
ions. A globally averaged water vapor concentration profile appropriate for 400 HLLV launches per year 
(see Sec. 6.1) was used in calculating the perturbed ionosphere, and the number densities of NO were 
fixed at ambient levels. HLLV- induced water vapor enhancements are appreciable (i.e., >20%) only in the 
upper D-region where water clustering to ions is inefficient. At lower altitudes (<80 ion) where clus- 
tering IS rapid, the computed water vapor concentration increase was <15%. This small increase in water 
vapor could not cause a marked change either in n^ or in the ionic composition, because there is 
already significant clustering at ambient water vapor levels. Although temperature changes alter the 
rates of cluster processes, the temperature changes caused by H 2 O (and NO) deposition are expected to 
have a negligible effect. 

The adopted ambient nitric oxide concentrations, [NO], and the computed enhancements, A[N0] due to NO 
eposited by reentering space vehicles are shown in Table 6. Two cases are shown: (1) one using A [NO] 

from the 1-D photochemical model and (2) one using similar data from the 2-D model. The limitations of 
the two models relative to NO predictions discussed earlier suggest that these two A [NO] are good esti- 
mates of upper and lower bounds on A[N0]. The computed values of A[N0] (the difference between col- 
umns 3 and 5, and column 1) are as large as 30% so that the predicted changes in electron concentration 
are more noticeable than the effects of water vapor alone. This is evident from Table 6 which indicates 
small de creases in electron concentration at an altitude of 60 km and above 85 km. These decreases are 
the result of increases in the dissociative recombination coefficients due to clustering of water vapor 
to ions. That decrease exceeds the increase in ionization rate due to enhanced NO. 


Table 6. Predicted Ionospheric Perturbations Due to Nitric Oxide 

Deposition 


Alt. . 
km 


100 

90 

85 

80 

75 

70 

65 

60 


Unperturbed Perturbed 

[NO], [e|, Case 1 Case 2 

cm“ cm"^ [NO] [e] [NO] [e] 


1.0 

00 

4.5 

(4) 

1.0 

(8) 

4.5 

(4) 

1.0 

(8) 

4.5 

(4) 

4.8 

(7) 

7.5 

(3) 

4.8 

(7) 

7.4 

(4) 

4.8 

(7) 

7.4 

(3) 

1.5 

(7) 

4.1 

(3) 

1.6 

(7) 

3.9 

(3) 

1.5 

(7) 

4.1 

(3) 

1.0 

(7) 

8.0 

(2) 

1.1 

(7) 

8.3 

(2) 

1.0 

(7) 

8.2 

(2) 

1.0 

(7) 

6.1 

(2) 

1.3 

(7) 

6.5 

(2) 

1.2 

(7) 

6.4 

(2) 

1.0 

(7) 

1.7 

(2) 

1.4 

(7) 

1.9 

(2) 

1.4 

(7) 

2.0 

(2) 

9.6 

(6) 

2.9 

(1) 

1.7 

(7) 

3.4 

(1) 

1.5 

(7) 

3.3 

(1) 

6.5 

(7) 

7.6 

(0) 

7.3 

(7) 

7.2 

(0) 

7.4 

(7) 

7.6 

(0) 


Collision 

frequency, 

a 


3.2 (4) 
1.9 (5) 
4.5 (5) 
1.1 ( 6 ) 
2.4 (6) 

5.3 (6) 

1.1 (7) 

2.2 (7) 


"^These values have been corrected for the energy dependence of the collision 
frequency (see Ref. 36). 

^ 1.0 ( 8 ) = 1.0 X 10 ®. 


The ionospheric D-region contributes significantly to the absorption of radio waves used in communica- 
tions. The absorption coefficient is given by (approximately) 


a = 5.31 X 10"® 


UV + (m ± cj^)2 


(52) 


where + and - signs refer to the ordinary and extraordinary modes, respectively, and where y, v, m, 
and n^ are the real part of the refractive index, electron-neutral collision frequency, angular fre- 
quency of the radio wave, and electron density, respectively (Ref. 34); wl is the product of the elec- 
tron gyrofrequency and the cosine of the angle between Earth’s magnetic field and the direction of the 
radio wave. When^ v >> w, a ^ n^/v; on the other hand, if v « oj, a - ngV. In either case, increased 
electron density in the D-region could lead to increased absorption of radio waves. To obtain a simple 
estimate of this effect, let us consider HF waves (3-30 MHz) for which the D-region absorption would be 
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given by 


, , 2.33 X 10-' 


n^v dh dB (mks units) 


(53) 


toward the upper portion of the range (Ref. 34) where the absorption is in dB and the integral is .over 
the entire D-region (assumed here to extend from 65 to 85 km). Using the values of n^ and v from 
Table 6, it is clearly seen that the D-region absorption of HF radio waves may be increased by about 10% 
during daytime. In comparison, much larger increases are encountered during large-scale geophysical 
disturbances (e.g., sudden ionospheric disturbances) which create serious communication problems. It is 
therefore doubtful whether the predicted 10% increase in D-region absorption would create a serious 
problem. Nevertheless, the question deserves more detailed examination. It should be pointed out that 
expression (53) does not apply at frequencies where m - v. Then a more complicated expression involving 
Dingle functions must be used (see Ref. 36). 


Propagation of extremely low frequency (ELF) waves is controlled by the low frequency, or dc, conductiv- 
ity of the air, particularly the conductivity scale height (Ref. 61). Our calculations show that the dc 
conductivity is hardly affected, the scale height even less so. Moreover, ELF communication links cover 
large horizontal distances, that is, latitudinal zones over which A [NO] and, consequently, changes in 
dc conductivity, may vary considerably in any case. The simple approach of Ref. 61 suggests a negligible 
effect of SPS on ELF propagation. 

We looked for but could find no evidence that metal atoms would be deposited by ablation during HLLV 
reentry. 

9,2 Short-Term Local Effects 


Initially the nitric oxide > produced by a reentering SPS vehicle is extremely localized along the entry 
trajectory. With the passage of time, this dense column of NO expands over wide regions. Table 7 pre- 
sents computed NO densities at various altitudes and at two times after the reentry, 6 hr and 1 day, 
obtained on the basis of the simplified NO dispersion model discussed in Sec. 2.3. The corresponding 
daytime equilibrium electron concentrations are also shown in these tables. Clearly, a significant local 
ionospheric absorption could occur near the reentry point. There is therefore a need to reexamine this 
possibility, using more accurate mathematical models. It is also necessary to consider micrometeorolog- 
ical factors, for example, local winds, which might force the initial column to evolve differently on 
different days. On the other hand, it seems unlikely that large regions would be affected by radio wave 
disturbances . 


Table 7. Nitric Oxide and Electron Concentrations Near Reentry 







Densities , 

, cm 

3 


km 


NO 



e“ 


Unperturbed 

6 

hr 

1 

day 

6 hr 

1 

day 

tiq from table 

100 

2.0 

(+9) 

6.0 

(+8) 

6.0 (+4) 

4.8 

(+4) 

4.5 (+4) 

90 

9.0 

(+10) 

9.0 

(+9) 

2.7 (+5) 

8.6 

(+4) 

7.5 (+3) 

80 

6.2 

(+10) 

1.2 

(+10) 

1.8 (+5) 

7.3 

(+4) 

8.0 (+2) 

70 

1.5 

(+11) 

2.5 

(+10) 

2.6 (+4) 

8.0 

(+3) 

1.7 (+2) 

60 

8.0 

(+10) 

1.5 

(+10) 

1.6 (+1) 

8.5 

(-K)) 

7.6 (+0) 


9. 3 Uncertainties 


Ionospheric D-region predictions are based on very uncertain physical systems. The uncertainties arise 
from a lack of information about the reactions occurring in the ionosphere, and of experimental values 
for certain parameters, such as rate constants and ionic mobilities. Observations of D-region electrons 
and ions often produce conflicting and inexplicable results. Thus, the accuracy of current ionospheric 
predictions is probably reasonable, but questionable. They are also limited by the accuracy of the pre- 
dictions of neutral atmospheric changes. The predictions for short-term ionospheric effects are expected 
to be more uncertain than the long-term predictions, because the factors controlling the spatio-temporal 
evolution of A [NO] are more complex. The nature of the predicted ionospheric effects, however, suggest 
that this problem may need to be Investigated more carefully. Our present results should therefore be 
regarded as preliminary estimates. 
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10. POSSIBLE CLIMATIC EFFECTS OF LARGE ROCKET OPERATIONS 

Some of the potential climatic effects of rocket clouds, for example, were mentioned in previous sections 
of this report. Here we discuss the more general mechanisms for climatic change. 

Changes in the mixing ratios of radiatively active minor constituents caused by HLLVs could conceivably 
alter atmospheric heating rates and produce climatic changes (climate here is defined narrowly as the 
time— averaged meteorological state at Earth’s surface). At least a 10— yr average would be required to 
obtain a representative climate that is not polluted by sampling error. Since compositional changes due 
to HLLVs are confined to levels above the midstratosphere, we must look for mechanisms that couple these 
higher levels with the troposphere. Such mechanisms can be radiative or dynamical. Radiative mechanisms 
involve perturbations of upper atmospheric composition or temperature that affect heating at the surface 
either through changes in the extinction of visible solar radiation or through changes in the transport 
irifrared radiation. Dynamical mechanisms involve high level circulation changes, including changes in 
reflection and absorption of planetary waves, that force changes in tropospheric winds, pressure, and 
temperature. 


Radiative mechanisms for climatic change have been examined in a number of radiative— convective one- 
dimensional models, including those of Ramanathan et al. (Ref. 62), Schneider and Coakley (Ref, 63), and 
Manabe and Wetherald (Ref. 64). Such models examine the globally averaged radiative balance for various 
atmospheric compositions, with dynamics essentially excluded, except for the constraint that vertical 
temperature gradients cannot exceed the adiabatic lapse rate. An examination of the results of these 
models indicates that the most important changes in composition associated with HLLVs, namely decreases 
in ozone and increases in H 2 O and NO^j would not have a significant effect on climate, at least as mea- 
sured by changes in surface temperature. For example, Ramanathan et al. (Ref. 62) computed a surface 
temperature increase of only 0.06°K for a 10% increase in stratospheric water vapor. Since HLLVs are 
expected to produce at most a 0.5% increase in water vapor, their globally averaged effect through this 
mechanism can be ignored. Ozone decreases should be even less important, since they are concentrated 
strongly at levels well below the maximum in ozone heating at 50 km. Even if ozone were to be decreased 
at all heights by the maximum value of 1%, surface temperature decreases would amount to no more than 
0.01°K (Ref. 62), This decrease would be made even smaller by the expected increase in NO^, which acts 
to warm the troposphere by absorbing solar radiation reflected from the surface. Thus, on a globally 
averaged basis, radiative mechanisms would be ineffective in causing significant climate changes asso- 
ciated with HLLVs, 

It is conceivable that the radiative effects could be enhanced by the formation of "corridors" of H 2 O 
and NOx enhancement (or ozone reduction). For example, if all the water vapor from the HLLVs were con- 
centrated in a 5® wide latitudinal band, water vapor increases of 10% within that band would be possible. 
However, this is not likely to result in a local 0.06®K increase in surface temperature, for horizontal 
heat transport will smear out any synoptic scale heating on a much more rapid time scale than the several 
months it takes to build up significant water vapor increases. In any case, it was shown in Sec. 6 that 
the establishment of a water vapor "corridor" enhancement is highly unlikely. 

It should be noted that cloud cover is an extremely important parameter in the assessment of climate — 
and a critically uncertain one as well. However, it appears that changes in cloud cover associated with 
HLLVs will be insignificant. In the mesosphere, noctilucent cloud incidence is not likely to increase 
greatly, as shown in Sec, 8 of this report. As for the troposphere, changes in water vapor will be less 
than 0.1%, hardly important for climatic change. Moreover, very little sulfur or particulate materials 
are expected from HLLVs, so that condensation or nucleation processes are not likely to be affected 
significantly. 

Dynamical interactions between the troposphere and upper atmosphere have traditionally been thought of 
as a "one-way street," that is, the troposphere performs work on the stratosphere, but does not respond 
^nrn to stratospheric motions. This is because stratospheric motions contain much less energy than 
tropospheric motions, which is due to the much lower mass of the stratosphere. The mean meridional 
circulation of the stratosphere, for example, would not have much effect on the troposphere since 
meridional wind speeds are comparable in magnitude to tropospheric wind speeds and thus much less ener- 
getic because of lower air density in the stratosphere. However, eddy motions, such as large-scale 
planetary waves, generally amplify as they propagate upward, maintaining a nearly constant energy per 
unit volume. If one allows for these waves to be reflected in the stratosphere by some mechanism, with 
the amplitude and phase of the reflected waves being governed by stratospheric conditions, it is pos- 
sible for the stratosphere to have a feedback effect on tropospheric planetary waves. In fact, as 
Charney and Drazin (Ref. 65) pointed out, the westerly jet of the winter upper stratosphere does act as 
just such a reflection mechanism. The implications of this for changes in climate are significant, for 
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if changes in the zonal jet could be reduced by compositional changes at high altitudes, planetary wave 
reflection characteristics would also be affected, thus affecting the phase or amplitude of planetary 
waves in the troposphere (Ref. 66). 

A number of investigations have explored this mechanism, including those of Avery (Ref. 67) and Geller 
and Alpert (Ref. 68). These studies focused on the linear planetary wave response of a number of zonally 
averaged atmospheres to topographic and diabatic forcing. Comparison of results for different cases 
could then show whether different types of changes in the stratospheric zonal wind had significant 
effects on tropospheric climate. The major result of the studies is that zonal wind changes of 20% or 
more must occur below about 35 km for planetary wave phase or amplitude changes at the surface to 
approach the Interannual variability (Ref. 68). Of course, this result does not take into account the 
highly nonlinear nature of stratospheric motions; still, it is a good indication of the magnitudes 
necessary to cause significant dynamically induced climate change. 

The implications of this result for SPS operations can be crudely assessed by making some estimates for 
the effects of HLLV water vapor emission. The radiative-convective model of Manabe and Wetherald 
(Ref. 64) suggests a 10®K decrease in stratospheric temperatures associated with a quintupling of water 
vapor . Scaling this down to an expected 10% maximum H 2 O increase if HLLV water vapor remains at low 
latitudes, we obtain temperature decreases of ~0.2®K. Comparing this with the equator-pole temperature 
gradient of about 40° (which is proportional to the zonal wind strength), it can be seen that the maxi- 
mum expected change in the zonal wind is only 0.5%, well below the threshold for significant effects on 
reflected planetary waves. As before, HLLV- induced ozone changes would have even less effect, since 
they are confined to levels above 40 km. 

Based on these estimates, we conclude that no significant effects on Earth* s climate should result from 
the projected launch schedule of HLLVs (400 yr”^). That is, we expect no temperature changes, locally 
or globally averaged, in excess of 0.1°K and no dynamical variations exceeding 5% of the interannual 
variability in planetary wave activity. 


1 1 . SUMMARY 

Based on detailed theoretical analysis, it is apparent that we have uncovered no serious environmental 
hazards associated with operating HLLVs in the middle atmosphere. On the contrary, we conclude that: 

1. "Corridor" effects, that is, enhancements of important constituents in a narrow latitude zone cen- 
tered on the launch point, are not likely to occur. The reason is that meridional transport is suffi- 
ciently rapid to spread the effluents fairly uniformly over the globe before a steady state is reached. 

2. With the possible exception of mesospheric nitric oxide, changes in neutral composition of the atmo- 
sphere below 90 km are not likely to be significant. Specifically, ozone decreases would probably not 
exceed 0.1% globally. 

3. Deposition of carbon dioxide in the stratosphere should not cause noticeable changes in the strato- 
spheric thermal balance; thus, the coupling to photochemistry is insignificant. 

4. The climatic effects of extended HLLV operations are probably negligible. The effects considered 
here include enhancements of stratospheric aerosol opacity due to SO 2 deposition by launch motors and 
alteration of planetary wave propagation characteristics in the stratosphere. 

5. Perturbation of HF, VLF, and ELF communication links by launch and reentry operations will be minimal 
compared to normal variability and natural disturbances of the ionosphere. The widespread ionospheric 
effects are small (VLF and ELF), and short-term effects are localized (HF) < 

On the other hand, there are some potentially important SPS effects that require further research: 

6. There may be a worldwide buildup of thermospheric hydrogen, perhaps a doubling of the hydrogen 
globally for continuous H 2 and H 2 O deposition above 70 km. The enhanced hydrogen in the magnetospheric 
ring current region may reduce the population of current-carrying protons through charge exchange. 

7. Extensive rocket contrails may form during HLLV launch and these may persist for a day or more, 
although extensive triggering of noctilucent clouds is highly unlikely. 
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8. The D-region electron number density and conductivity may be greatly enhanced locally during atmo- 
spheric reentry. However, the modified region must be large if VLF and ELF propagation were to be 
affected. 


On that note we conclude this report, hoping that future work will resolve the questions raised in 
items (6) and (8). 
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